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Abstract: This paper introduces in brief the traction system of a permanent magnet 
electrodynamic suspension (EDS) train. The synchronous traction mode based on long 
stators and track cable is described. A speed and position detection system is 
recommended. It is installed on board and is used as the feedback end. Restricted by the 
maglev train's structure, the permanent magnet electrodynamic suspension (EDS) train 
uses the non-contact method to detect its position. Because of the shake and the track 
joints, the position signal sent by the position sensor is always aberrant and noisy. To solve 
this problem, a linear discrete track-differentiator filtering algorithm is proposed. The 
filtering characters of the track-differentiator (TD) and track-differentiator group are analyzed. 
The four series of TD are used in the signal processing unit. The result shows that the 
track-differentiator could have a good effect and make the traction system run normally. 

Keywords: permanent magnet electrodynamic suspension train; position sensor; speed and 
position detection; track-differentiator; filtering character 



1. Introduction 

A permanent magnet electrodynamic suspension (EDS) train uses a linear motion actuator to 
achieve traction [1]. Its current transformer is placed on the ground. The long (primary) stators of the 
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linear synchronous motor are laid on the middle of the tracks. The permanent magnet for traction is 
installed at the bottom of the vehicle which faces the long stators. This is shown in Figure 1. This kind 
of traction mode has high power factor and high efficiency the characteristics, allowing it to attain high 
acceleration and speed. The coils of the long stators can produce a travelling wave magnetic field. To 
make the traveling magnetic field synchronize with the magnetic field produced by the permanent 
magnet, the traction system needs to detect the position of the secondary, only so the current can be 
modulated to control the output force. Then the detection of the secondary position is extremely 
important. The error detection will not only make the vehicle run abnormally, but also cause serious 
damage to the correlated equipments [2-4]. 



Figure 1. The long stator and track structure. 




The synchronous traction system of a permanent magnet EDS train comprises the speed and 
position detection system, the radio unit, ground traction system and traction power module. It is 
shown in the Figure 2. In the synchronous traction system, there are two methods to detect position of 
the motor's secondary. One method is to detect the counter electromotive force. The other is achieved 
by the position sensors. The former one is called sensorless. When the train runs at a very low speed, 
the sensorless method cannot work very well, because of the weak counter electromotive force and 
strong noise. During the slow running period, the position of the train is determined by sensors. This 
kind of method is also used in electric rotating motors. The electric rotating motor uses a photoelectric 
encoder, rotary transformer or inductive rotation sensor to detect position or speed. The character of 
the inductive rotation sensor is that the distance between detection probe and speed gear is non-contact 
and unchanged. But the same usage in a permanent magnet EDS train has its own special problems. 
The position sensor is a kind of on-board equipment. Its detection gap fluctuates with the levitation 
gap. In addition, its faced stack long stators do not have a continuous conductor surface, because the 
long track stators are spliced with joints of different sizes. Also the sensors are easily affected by all 
kinds of magnetic field, including the traveling magnetic field and the suspension magnetic field. 
These problems are different from the speed detection of the electric rotating motor, so this paper 
examines the sensors' signal processing, dealing especially with error signal. 
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Figure 2. Composition of the synchronous traction system, 
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2. Principle of Speed and Position Measurement Based on the Long Stators 

For a maglev train, the speed and position detection system is necessary and crucial [5]. It carries 
out the task of detecting the speed and position of the train. The operation control system also takes 
advantage of it to achieve centralized control. In order to make the speed and position detection system 
reliable and safe, a convenient structure is needed [6,7]. This is shown as Figure 3. 

Figure 3. Structure of the speed and position detection system. 
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The position sensor's detection object is the long stator which is used as the train's traction. The 
long stator is composed of silicon steel sheets. It has the alveolar structure shown in Figure 3. The 
detection coils of the sensor face the alveolar structure. When a sensor moves along the long stators, it 
can distinguish the teeth and the slots. By counting the passed alveolar structure, the train's position is 
detected [8,9]. Generally, the position sensor has the three functions: 

(1) Alveolar counting: The position sensor counts the passed alveolar structure, and then sends the 
result to the signal processing unit. 

(2) Phase detection: The detection accuracy can only reach one alveolar structure's length which 
can hardly fulfill the requirements of the traction system, so the position sensor needs to give a 
subdivision. The phase signal between zero degrees to sixty degrees represents the length of 
one alveolar structure. 

(3) Judging direction: The direction signal is not only the foundation of the sensor's alveolar 
counting, but also the precondition of the unit of signal processing to judge the train's direction. 
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The signal processing unit is the core component of the speed and position detection system. Each 
system has one signal processing unit, which is used to receive and handle the sensors' data. The signal 
processing unit sends the data processing results to the traction system in terms of the agreement. To 
some extent, it is the bridge between the traction system and all kinds of sensors. Generally, the signal 
processing unit has the following main functions: 

(1) Receiving the sensors' position signal, speed signal and direction signal; 

(2) Data filtering and information fusion in accordance with the requirements of the traction system; 

(3) Communicating with the protocol processing unit. The processed position and direction data is 
sent to the traction system by the protocol processing unit. 

The data that the signal processing unit sends to the traction system includes the following content: 

(1) Magnetic pole phase (zero degree to three hundred and sixty degrees): The magnetic pole phase 
signal of is used for the traction, where every sixty degrees represents one alveolar structure. It 
is calculated by the alveolar counting signal and phase detection. 

(2) Magnetic pole counting (0 to 65,535): As shown in the Figure 4, each time the train passes six 
alveolar structures, the magnetic pole counting signal increases by one or decreases by one in 
terms of the direction. 



Figure 4. The position sensor and the long stator. 
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3. Characteristics of the Discrete Tracking-Differentiator Filter 

3.1. The Linear Discrete Track-Differentiator 



The classical differentiator makes use of an inertia element to track the input signal with a time 
delay. It has the following form [5]: 

U{t)-U{t-T) 



u\t) ^ - 



(1) 



But the classical differentiator also amplifies the noise, so another approximate form is always used: 

(2) 



, u(t -T,) -u(t -T.) ^ 
u\t) - ^ ^,0 <T,<T^ 



where, the signals u{t-T^) and u{t-T^) are obtained from the inertia element. Then, the noise 
amplificatory effect can be reduced. Then there is the transfer function of the differentiator: 

1 



:^(-^ 



)u 



(3) 
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If the time-delay parameters and are equal to co , the formula (3) can be turned into the 
following expression: 



0) s 



Y{s) 



(4) 



(5) 



It can be discretized into a two order equation of state: 

{x^{k + \) = x^{k) + h-x^{k) 

\x^{k + \) = (k) + h' i^-o/ (x^ (k) - u(k)) - l^cox^ (A:)) 

The formula (5) is the discrete form of the linear track-differentiator. The variable x^ (k + 1) tracks 
the input signal u(k) .The parameter ^ is the damping factor and co is the quickness factor. The bigger 
the dereferencing of co is, the better x^(k + l) approximates u(k) . And then x^ik + l) is the 
approximate value of u(k) . However, the linear track-differentiator's noise immunity is very limited. 
In order to suppress the noise and track the input signal quickly, this paper adopts a new nonlinear 
track-differentiator filter to deal with the input signal. The track-differentiator this paper uses is 
different from Han's one [10]. The retum-to-zero strategy near the switching curve is proposed [11]. 

3.2. Signal Filter Based on the Nonlinear Track-Differentiator 

First of all, we consider the problem of a second-order continuous-time system: 



X, 




0 1 


X, 
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^2 




0 0 


^2 
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(6) 



The initial state is X(^o) = the permissible control domain is \u\ < r .The problem is that how 

to find the time optimal control to make the system performance index j = [ dt minimum. 

There exists one optimal control strategy. It can make arbitrary point in the phase plane reach the 
origin using at the most once switching point and can spend the shortest time in doing so. The strategy 
is called bang-bang control. All of the control switching points compose the switching curve: 



( 1 ' 2 ' ) 1 2 I 2 



(7) 



According to the optimal control theory, the optimal solution of the optimal control problem above 
is as follows: 



u(t) = ' 



■ sign(r[x^,X2,r)^ ,T[x^,X2,r) ^ 0 



I r •^/g^(x2) , r[x^,X2,r) = 0 , ^2 ^0 

If (t) is substituted by (t) - v{t) , then the switching curve is as follows: 

r{x„X2,r) = x,(t)- v(t) + -^^2 (0 • |x2 (0| 

The continuous-time track-differentiator is obtained: 



(8) 



(9) 



-r • sign 
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If the arbitrary point Mix^^^x^^) is not in the switching curve, is used to represent the time 
thatM(xio,x2o) spends to reach the switching curve, is calculated as follows: 



r 2 



(11) 



If the arbitrary point M(x^q,x2q) is in the switching curve, is used to represent the time that 
M(xio,x2o) spends to reach the origin, is calculated as follows: 



t. = 



(12) 



In order to get the fastest control function, the step-size unchanged strategy is adopted. The 
paper [6] proposes the discrete form of nonlinear track-differentiator: 



u(k) 

b 



= fast(x^(k)-v(k),c^X2(k),r,CQh) 



x^ {k + \) = x^ {k) + h'X2 {k) + r~h^ 'u{k)-^b'h^ 
X2{k-\-l) = X2{k)-\-r'h'u{k)-^b'h\k = 0X2"'. 



(13) 



where x^(0) = x^^, X2(0) = X20 . 

The termr is the quickness factor. It decides the tracking time, h is the sampling period. A suitable 
h makes the filtering effect better. If h is to too large, the phase of the follow-up signal may suffer 
from errors. Unless stated otherwise, this paper chooses one millisecond as sampling period, q is the 
damping factor and is the filtering factor. The nonlinear function fast(x^,x2,r,h) is inferred by the 
rules as follows: 

(1) When M (x^(k),X2(k)) hasn't reached the switching curve, meanwhile t^> h , it has the 
following expression: 

{b = 0 

(14) 



u(k) = -sign Xj (k) H X2 (k) • 1x2 (k)\ 

y 2-r 

(2) When M {x^{k),x^{k)) hasn't reached the switching curve, meanwhile t^<h , it has the 
following expression: 

{b = Q 

^ ^ (15) 



1 X2(A:)-^ 1 
u(k) = — + — + - 

2 r-h 2 



1 + 



4-X2(^)_^8-Xi(^) 



r-h 



r-h' 



where s = sign(r[x^ , X2 , ^)) = signi Xj H — — X2 -1x2 

V 2-r 

(3) When M (x^(k),x2(k)) hasn't been in the switching curve, meanwhile > A , it has the 
following expression: 



b = 0 

u(k) = -sign[x2(k)) 



(16) 



(4) When M (x^(k),x2(k)) hasn't been in the switching curve, meanwhile t^<h ^ it has the 
following expression: 
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TP (17) 

u(k) = a + bh 

where a- H^'-2^k)M^x,{k)) 

r-h^ 

In general, the function of fast has the form as this: fast{x^-v{k\x^,r,CQh) , which just needs some 
substitution work for the four rules above. 

On the one hand, this track-differentiator can quickly track an input signal and can produce a good 
differential signal. However the differential signal is not used in this paper. On the other hand, the 
output signal's amplitude has almost no attenuation and its time delay is small. Because of the little 
attenuation of the signal amplitude and linear relations between the time delay and signal frequency, 
the true value can be separated effectively from signals mixed with noise, so this track-differentiator 
has the capability to act as an effective filter [11-13]. The relevant proof can been found in [13]. 

3.3. Signal Filter Based on Track-Differentiator Group 

As far as the single track-differentiator, the output signal k is a kind of function on the filtering 
factor Co, damping factor ci, quickness factor r and time-delay factor r. The output signal of 
track-differentiator has a time delay. The solution is to modulate the correlated parameters, but 
tracking errors will then increase. The track-differentiator group is a kind of method to relieve the 
contradiction error and time delay. It is shown as Figure 5. The track-differentiator filter is called TD 
for short. 

Figure 5. Comparison between track-differentiator and track-differentiator group. 
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where V{t) is the signal to be tracked. The tracking error is marked as ^{t). It is caused by the filtering 
system and the input signal. Considering the tracking error and time delay, the output of TD is as follows: 

V(t) = V(t-r) + ^(t) (18) 

With regard to the system with one TD, Taylor expansion oi V{t- t) about t is obtained: 

= K„,r(0(-,).a££(_,,. ,a£zpi):(_.). (19) 

In the case of a TD of series length m, the output of TD. is the input of TD.^^ , if necessary. The 
input signal of the track-differentiator group passes by a number of TD, an output with a nice filtering 
effect will be obtained. So the track-differentiator group may overcome the drawback of one TD's 
weak filtering. Meanwhile the tracking error is considered to approach to zero. But the phase 
compensation comes to be the principal problem. One TD may bring the time delay of T , then a 
length of m series of TD may result in m times of r . Some sort of phase compensation is needed. 
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Firstly, the track-differentiator group's input-output model is given as follows: 

\V,{t) = V{t-T) 

V,{t) = V{t-2t) 
V,{t)^V{t-4T) 
y„{t)^V{t-mf) 

Like in formula (19), Taylor expansion n of F(? - r) about t is obtained: 

V"(t\ V''(t\ 
V(t- iT) = F(0 - ViO-iT + -^'{irf ^{iry + R{t, r) 

If m = 4 , the high order term is omitted, it has the following form: 



1 -2r 
1 -3r 
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(20) 
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(22) 



To make it further, the parameter of time delay T is drawn out from the matrix of coefficients, and 
then the matrix equation has the following form: 



1 -1 1 -i 
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7. 



To get expression which is relative every TD' output, the matrix of coefficients is inverted as follows: 

(24) 

From the matrix equation of (24), the approximation of input signal V(t) is obtained and it also 
serves as the output of track-differentiator group: 

V(t) = 4V, (t) - 6V, (0 + 4^3 (0 - V, (t) (25) 

Considering the generalized case, the tracking signal with phase compensation is as follows: 

m 

K=T.(^miK (26) 

When the variable m is assigned different value, the coefficient of a^. makes a difference. 
Theoretically, as the value of m increases, the calculations needed will be considerable. In actual use. 
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the order of the track-differentiator should be chosen based on the system's requirement and 
computational capabilities. 



Table 1. Coefficients a, 
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4. Signal Process for System of Position and Speed Detection 

4.1. Description of the Sensor's Abnormal Signal 

In theory, the signal of phase detection increases or decreases in terms of the train's running 
direction between zero degrees to sixty degrees. In the forward running condition, as the phase 
detection signal changes from sixty degrees to zero, the alveolar counting signal increases by one. In 
the backward condition, as the phase detection signal changes from zero degree to sixty, the alveolar 
counting signal decreases by one. But in the real circuit processing, the sine wave obtained by 
synchronous demodulation is sent to a different post-processing circuit, and then the change of phase 
signal and alveolar counting signal are not synchronized. If the signal processing unit simply integrates 
them, there will be a mutation as shown in Figure 6. The amplitude is about sixty degrees. 

Figure 6. The magnetic pole phase because of asynchronization. 
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As the position sensor moves along the long stators, the inductance of its coil will change. In short, 
the position sensor can detect the change of different metal's inductance. Meanwhile the position 
sensor is installed at the bottom of the train, so its detection gap is changeable with the suspension 
height. Then the fluctuating gap may also make the inductance change, which makes the vertical and 
horizontal changes coupled. In other words, the sensor can hardly distinguish between vertical and 
horizontal changes, so the linearity of its position detection will be affected. 
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For prevention of thermal expansion and contraction and also for convenience when laying out the 
long stators, the long stator track is composed of numerous stators, so there are many joints. Usually, 
the length of the joint is seventy to eighty millimeters. This is shown in Figure 7. There is nothing 
except track cable in the joint. The existence of joints breaks the continuity and the periodicity of track 
inductance and increases the difficulty of detection. As the position sensor is passing the track joint, 
the position signal will be more or less aberrant. Figure 8 shows the affected phase. 

Figure 7. The joint of long stators. 
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Figure 8. The aberrant signal affected by joint. 
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4.2. Pre-Treatment of the Position Sensors 



The signal processing unit receives the phase and alveolar counting signal to locate the train. Then 
the train's position is equal to the phase signal plus the alveolar counting signal. To avoid the low-pass 
effect of the track-differentiator, the signal processing unit needs to pre-treat the signals. 

The frequency spectrum of the position signal must be within the track-differentiator's frequency 
band, so the high-frequency noise can be filtered. Considering the relation of phase and alveolar 
counting, it can be easily found that the phase is a subdivision of alveolar counting. The phase signal 
can be treated as the low bit while the alveolar counting is the high bit, so the saw-tooth wave is turned 
into a ramp wave. Then they can be integrated as follows: 

Firstly, considering the range of phase signal x^^ , it can be normalized into the range of zero to ten. 
The method is as follows: 



/;; 



60 



(27) 
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Secondly, the normalized result p can be integrated with the signal of alveolar counting;;. It is 
marked as formula (28). Especially under uniform velocity conditions, the integrated result will be 
similar to a ramp wave. Figure 9 shows the forward phase signal and the integrated result. Figure 10 shows 
the backward phase signal and the integrated result. Note that the integrated signal is dimensionless. 

P = x,,+y\Q (28) 



Figure 9. The phase detection and integrated signal in the forward direction. 
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Figure 10. The phase detection and integrated signal in the backward direction. 
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4.3. Application of the Track-Differentiator in Speed and Position Detection System 



The sample period of the discrete track-differentiator is one millisecond, and the quickness factor 
r is 1,500. Figures 11 and 12 show the comparison of the input signal and the output signal. It can 
be seen that the track-differentiator can effectively smooth the input signal and remove the impulse 
interference. 

From comparison of the input curve and the following curve in Figure 12, the filtering results in the 
time delay just as the analysis of Section 3.3. Even though the filter produces a nice smooth effect, the 
output signal cannot be used because of the time delay. 
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Figure 11. Comparison of the input and the output signal. 
5250 
5200 
5150 
5100 
5050 
5000 
4950 

2.7 2.75 2.8 2.85 

Tinne(nnillisecond) ^ 

Figure 12. Partial enlarged drawing of Figure 11. 
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Meanwhile the damping factor q can be decreased to make time delay smaller. But if q gets 
smaller, the filtering effect will get worse. It is shown in Figure 13 that furthermore the quickness 
factor r and the filtering factor are modulated, the conflict of filtering and time delay still cannot 
be resolved. 
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4.4. Application of the Track-Differentiator Group in Speed and Position Detection System 

Considering the same noisy conditions of Section 4.3, four series of TD are used. Each TD has the 
same quickness factor, filter factor and damping factor, where r = 1500,CQ=2,q=1.5. The output of 
the track-differentiator group is shown in Figures 14 and 15. It can be seen that the output signal can 
follow the input signal nearly without time delay. Also Figure 15 shows that the pulse noise has been 
eliminated and the aberrant signal has been repaired. 

Figure 14. Comparison of the input signal and the output signal. 
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Figure 15. Partial enlarged drawing of Figure 14. 
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According to the introduction of Section 2, the signal processing unit needs to send the magnetic 
pole phase to the traction system, so it is recommended to recover the integrated signal to the phase 
signal and alveolar counting signal. Because the recovered phase signal is processed by the 
track-differentiator group, it is smoother and has less noise. This is shown in Figure 16. 
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Figure 16. Comparison of the phase detection and processed phase signal. 
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5. Conclusions 

This paper analyses the reason that the phase detection signal and alveolar counting signal of the 
position sensor are unsynchronized. Aiming at eliminating the noise and overcoming the effect of the long 
stator joints, a new non-linear discrete track-differentiator is proposed. Meanwhile the track-differentiator 
group is also designed and analyzed. Especially, a four series TD is used in the speed and position 
detection system systems of a permanent magnet electric suspension train. In the train's actual 
operation, it has played a good effect. 
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